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ABSTRACT 

Upcoming or future deep galaxy samples with wide sky coverage can provide independent mea¬ 
surement of the kinematic dipole - our motion relative to the rest frame defined by the large-scale 
structure. Such a measurement would present an important test of the standard cosmological model, 
as the standard model predicts the galaxy measurement should precisely agree with the existing pre¬ 
cise measurements made using the CMB. However, the required statistical precision to measure the 
kinematic dipole typically makes the measurement susceptible to bias from the presence of the local- 
structure-induced dipole contamination. In order to minimize the latter, a sufficiently deep survey 
is required. We forecast both the statistical error and the systematic bias in the kinematic dipole 
measurements. We find that a survey covering ~ 75% of the sky in both hemispheres and having ^ 30 
million galaxies can detect the kinematic dipole at 5 ct, while its median redshift should be at least 
Zmed ~ 0.75 for negligible bias from the local structure. 


1. INTRODUCTION 

Measurements of the motion of our Solar System 
through the cosmic microwave background (CMB) rest 
frame represent one of the early successes of precision 
cosmology. This so-called kinematic dipole corresponds 
to a velocity of (369±0.9) kms“^ i n the direction (Z, b) = 
(263.99° ± 0.14°, 48.26° ± 0.03°) (iHinshaw et al.l[200l . 
The kinematic dipole has even been detected (though 
not as precisely measured) by observing the relativistic 
aberration in the CMB anisotropy that it causes, which 
is detec ted via the coupling of high CMB mnltipoles in 
Planck (lAghanim et al.lr20i^ . 

Independently, the past few decades have seen signifi¬ 
cant progress in measuring the dipole in the distribution 
of extragalactic sources. The contribution of our motion 
through the large-scale structure (LSS) rest frame - the 
kinematic dipole - also leads to relativistic aberration, 
this time of galaxies or other observed LSS sources. We 
define the dipole amplitude via the amount of its “bunch¬ 
ing up” of galaxies in the direction of the dipole 
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pole £ = 1 in all extant LSS survey s. It has been mea¬ 
sured and reported e ither explicitly (iBaleisis et al.lll998l 
Blake fc Walll [200^ iHiratal 120091: iGibelvou fc Huterer 


2012 : iFernandez-Cobos et alJ 20141: iRubart fc Schwarz 

2013 : _ Yoon et alJ I20l4 lAppleby fc Shafielool 12014 

Alonso et al.ll201^ . or as part of the angular power spec¬ 


trum measurements. No LSS survey completed to date 
therefore had a chance to separate the small kinematic 
signal from the larger local-structure dipole contamina¬ 
tion due to insufficient depth and sky coverage. This will 
change drastically with the new generation of wide, deep 
surveys. 

Standard theory based on the adiabatic initial pertur¬ 
bations predicts that the kinematic dipole measured by 
the LSS should agree with the one measured by the CMB. 
Detection of an anomalously large (or small) dipole or the 
disagreement of its direction from that of the CMB dipole 
could indicate new physics: for example, the presence of 
superhorizo n fluctuations in the presence of isocurvature 
fluctuations (iTurneil Zibin fc Scottll200Mlltoh et ahl 


l20irt lErickcek et al.l 12008 ). Clearly, a kinematic dipole 
detection and measurement represent an important and 
fundamental consistency test of the standard cosmologi¬ 
cal model. 


where N is the galaxy number in an arbitrary direction 
n, d is the dipole direction, and e is random noise. The 
dipole amplitude A is approximately (but not exactly) 
equal to our velocity through the LSS rest frame in units 
of the speed of light; the precise relation is given in the 
following section. 

However, the dominant contribution to the LSS dipole 
is typically not our motion through the LSS rest frame, 
but rather the fluctuations in structure due to the finite 
depth of the survey. The dipole component of the lat¬ 
ter - the s o-called “local-structure dip ole” in the nomen¬ 
clature of IGibelvou fc Hute^ (|20I2ll - has amplitude 
A ^ 0.1 for shallow surveys extending to Zmax 0.1, 
but is significantly smaller for deeper surveys. The 
local-structure dipole is the dominant signal at multi- 


2. METHODOLOGY 
2.1. Theoretical signal 

The expe cted LSS kinematic dipole signal amplitud e 
is given by (jBurles fc RaDDaDortll2006t lltoh et aLll201fll l 

A = 2/3 = 2[1-L1.25x(1-p)]/3 (2) 

where (3 = v/c = 0.00123 (assuming the CMB dipole). 
The contribution 2/3 comes from relativistic aberration, 
while the correction [1 -|- 1.25a;(l —p)] corresponds to the 
Doppler effect; here x is the faint-end slope of the source 
counts, X = d\ogiQ[n{m < miini)]/dTOiim, and p is the 
logarithmic slope of the intrinsic flux density power-law, 
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Clearly, the parameters x and p depend on the pop¬ 
ulation of sources selected by the survey, and on any 
population drifts as a function of magnitude. We now 
estimate these parameters - note also that we only need 
the quantity A to set our fiducial model, so very precise 
values of th e population parameter s are not crucial for 
this paper. iMarchesini et al.l (|2012D find that the faint 
end of the V-band galaxy luminosity function does not 
vary much over the redshift range 0.4 < z < 4 and is 
equal to, in our notation, x = 0.11 ± 0.02. Moreover, 
for optical sources the flux density slope p varies signifi¬ 
cantly with the age of the source, but in the infrared it is 
more consistent, with mea s urements indica ting p ^ 0 
(|White fc Maiumdai] 12004 iMo et all I2010D . Here we 
adopt p = 0. Applying all these values to Eq. © , we 
get 

A ~ 0.0028 (expectation from CMB). (3) 

While the actual value of the kinematic dipole is of course 
unknown prior to the measurement, standard cosmology 
theory predicts it takes this value, plus or minus 0(20%) 
changes depending on the source population selected. We 
adopt Eq. ([3]) as the fiducial amplitude. 

The fiducial direction we adopt is the one of the best- 
fit CMB dipole, {l,h) = (263.99°, 48.26°). Note, however, 
that the results may vary depending on the relative ori¬ 
entation between the actual dipole direction and the cov¬ 
erage of the observed sky. Finally, note that bias (of the 
galaxy clustering relative to the dark matter field) enters 
into the contamination of the kinematic dipole measure¬ 
ments, but not the signal. The former quantity - the 
local-structure dipole - is linearly proportional to the 
bias b. Therefore, the bigger the bias, the more contam¬ 
ination the local-structure dipole provides for measure¬ 
ments of the kinematic effect. In this work we assume 
bias of 6 = 1. Note that the kinematic signal itself, be¬ 
ing due to our velocity through the LSS rest frame, is 
independent of bias. 


2.2. Statistical error 


Rewriting Eq. m somewhat, the modulation in the 
number of sources is given at each direction ri can be 
written as 
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where x = {dx,dy,dz,ki, ...,kM) is the vector of the 
three dipole component coefficients in the three spa¬ 
tial coordinates, plus coefficients corresponding to other 
multipoles (one for the monopole, five for that many 
components of the quadrupole, etc), as well as any 
desired systematic templates. The vector T(n) = 
(ria,, riy, n^, ti(n),..., tM(?T-)), contains the three dipole 
unit vectors (with nx+ny + nl = 1), plus M additional 
spatial patters for all templates included. Note that the 
choice of the fiducial values of the non-dipole template 
coefficients ki is arbitrary, since we will fully marginalize 
over each of these, effectively allowing ki to vary from 
zero to plus inf inity. The op timal estimate of x is given 
by X = F~^g (|Hiratal [200911 . where the components of 
the vector g are gi = f Ti{h)5N^{n)d^h and the best- 
fit dipole dbest is given by the first three elements of x. 


Here the Fisher matrix F is given by 

F,j =N^ Jt, [n)Tj {h)d^n, (5) 

where = dN / dH, is the number of galaxies per stera- 
dian and H is a solid angle. Note that the Fisher in¬ 
formation is proportional to the number of sources, and 
unrelated to the depth of the survey. It is therefore the 
number of sources, together with the sky cut (not just 
the fraction of the sky observed /sky but also the shape of 
the observed region relative to the multipoles that need 
to be extracted) that fully determines the statistical error 
in the various templates including the dipole. 

The Fisher matrix contains information about how well 
the three Cartesian dipole components, as well as the 
multipole moments of all other components, can be mea¬ 
sured in a given survey. Our parameter space has a 
total of Mpar = (^max + I)^ parameters, where imax is 
the maximum multipole included to generate the tem¬ 
plates (see below for more on the choice of ^max)- With 
this Fisher matrix in hand, we then marginalize over the 
M = Mpar — 3 non-dipole parameters, using standard 
Fisher techniques, to get the 3x3 Fisher matrix describ¬ 
ing the final inverse covariance matrix for the dipole com¬ 
ponents. Finally, we perform a basis change, converting 
from Cartesian coordinates {dx-,dy,dz} to spherical co¬ 
ordinates {A, 9, (j)} (where A is the amplitude of dipole), 
by using a Jacobian transformation to obtain the desired 
3x3 Fisher matrix in the latter space, 4^^® fore¬ 

casted error on A is then given in terms of this matrix 
as 

a{A) = (6) 

In a realistic survey with partial-sky coverage, the pres¬ 
ence of other multipoles (monopole, quadrupole, etc) will 
be degenerate with the dipole, degrading the accuracy in 
determining the latter. We have extensively tested for 
this degradation, in particular with respect to how many 
multipoles need to be kept - that is, what value of ^max 
(and therefore M) to adopt. We explicitly found that the 
prior information on the “nuisance” (7^, corresponding to 
how well they can be (and are being) independently mea¬ 
sured, is of key value: once the prior information on the 
Ci - corresponding to cosmic variance plus measurement 
error - is added, very high multipoles are not degenerate 
with the dipole. Our tests show that keeping all multi¬ 
poles out to irnax = 10 is Sufficient for the dipole error to 
fully converge. We also experimented with adding addi¬ 
tional individual templates ti (n) corresponding to actual 
sky systematics and with modified coverage (correspond¬ 
ing to e.g. dust mask around the Galactic plane), but 
found that these lead to negligible changes in the results. 
Moreover, we envisage a situation where the maps have 
already been largely cleaned of stars by the judicious 
choice of color cuts prior to the dipole search analysis. 
For these two reasons, we choose not to include any ad¬ 
ditional systematic templates in the analysis. 

2.3. Systematic bias 

The local-structure dipole diocai will also provide a con¬ 
tribution to the kinematic signal dkin that we seek to 
measure. The observed dipole in any survey will be the 
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error ellipsoid 

Vkin 



Fig. 1. — Sketch of the problem at hand: we would like to mea¬ 
sure the kinematic dipole i/kin, whose error (represented by a cyan 
ellipse) can be calculated given the number of extragalactic objects 
and the sky coverage. The LSS local dipole, Floral i provides a bias 
in this measurement. For a survey deep enough (and depending 
somewhat on the direction of its uiocaOi bias in the measurement 
of ffkin will be smaller than the statistical error. 

sum of the two contributions; 

dobs — dkin ~1“ diocal- (7) 

Without any loss of generality, we consider the kinematic 
dipole as the fiducial signal in the map, whose errors are 
therefore given by the Fisher matrix worked out above. 
We consider the local-structure dipole to represent the 
contaminant whose magnitude, ideally, should be such 
that the resulting observed dipole dobs is still within the 
error ellipsoid around the kinematic dipole direction and 
amplitude. This is illustrated in Fig. 1. 

We now quantify the systematic bias, due to the lo¬ 
cal structure, relative to statistical error in the measure¬ 
ments of the kinematic dipole. First note that we are in 
possession of the (statistical) inverse covariance matrix 
for measurements of the kinematic dipole, which 

is already fully marginalized over other templates. The 
quantity 

Ax^(diocal) = (dobs - dkin)^-F(““®)(dobs “ dkin) 

= dLal^^(3x’?)dlocal (8) 

then represents “(bias/error)^” in the kinematic dipole 
measurement due to the presence of the local-structure 
contamination. This chi squared depends quadratically 
on the expected local-structure dipole, and is therefore 
expected to sharply drop with deeper surveys which have 
a lower |diocai|, as we find in the next section. With three 
parameters, requiring 68% confidence level departure im¬ 
plies = 3.5. We therefore require that, for a given 
survey, the local-structure dipole magnitude and direc¬ 
tion are such that the value in Eq. (|5]) is smaller than 
this value. 

3. RESULTS 

For a fixed sky cut and number of sources in the survey, 
we first calculate the error on the amplitude of the dipole 
a{A). In the left panel of Fig. [5] we show errors as a func¬ 
tion of the number of galaxies in a survey. As previously 
noted, this statistical error does not depend on the depth 
of the survey, but does depend on both /sky and the shape 
of the sky coverage. Here we show results for an isolat¬ 
itude cut around the equator of ±15 deg, and ±30 deg 
and isolatitude cap-shaped cuts of 90 deg (i.e. half the 
sky removed) and 60 deg (i.e. leaving out a circular re¬ 
gion around a pole). Note that the Galactic ±15deg cut 
and the cap cut of 60 deg both have /sky = 0.75, while 


the Galactic ±30 deg and the cap 90 deg cuts both have 
/sky = 0.5. The results will also depend on the fiducial 
amplitude of the dipole, and here and throughout we 
assume the CMB-predicted value of A = 0.0028. Even 
for a fixed /sky of the survey, the cut geometry clearly 
matters, and the Galactic-cut cases have a smaller er¬ 
ror in the dipole amplitude due to symmetrical covering 
of the two hemispheres. For Galactic ±15 deg case, i-a 
and 5 -ct detections are easily achievable, requiring only 
A^gai = 9 X 10® and 3 x lO”^ objects, respectively. Note 
that if the actual direction of the LSS kinematic dipole 
deviates from the assumed dipole direction (the CMB 
direction), the result changes. For a 5-sigma detection 
and the same ±15 deg isolatitude cut, a dipole pointing 
along toward a Galactic pole, which is the best-case sce¬ 
nario, only requires 8 million objects; if instead the dipole 
points toward the Galactic plane, then 70 million objects 
are required. For the Galactic ±30 deg cut, the 3-cr de¬ 
tection is more challenging since it requires having over 
A^gai = 10® sources. The cap 60 deg cut mostly follows 
the trend of the Galactic ±30 case. Lastly, the 90 deg 
cap cut cannot detect the signal even at the 1-sigma level 
and with Agai = 10®. We conclude that dual-hemisphere 
sky coverage is crucial in the ability of the survey - or a 
combined collection of surveys - to detect the kinematic 
dipole. 

The right panel of Fig. [2] shows the systematic bias in 
the dipole measurement due to the presence of the local- 
structure contamination, showing the quantity defined in 
Eq. ([8]) . Because diocai has an a-priori unknown direction 
and its amplitude changes according to the depth of the 
survey, the systematic error is a function of direction of 
diocai and the depth of the survey. Therefore, we choose 
to plot Ay® averaged over all directions of diocai- 

To calculate the amplitude of diocai, we 
model the radial distribu tion of obj ects as 
n{z) = z®/(2zq) exp(—z/^o) (|Hutererl 1200^ . where 
the parameter zq is related to the median redshift as 
Zo = Zmed/2.674. A deeper survey (larger Zq) has a 
smaller local-structure dipole. Note that one could 
additionally cut out low-redshift objects in order to 
further reduce the contamination from the nearby 
structures, as well as the star-galaxy confusion; we have 
not done that here. 

The dashed lines in the right panel of Eig. [2] represent 
the cases when Zmed = 0.75 and the solid ones are when 
-Zmed = 1-0. Since Ay® is inversely proportional to the 
statistical error squared, the best cases in the left panel 
of Eig. [2] have a larger bias in the right panel. In par¬ 
ticular, the more galaxies the survey has, the more it 
is susceptible to systematic bias (for a fixed depth and 
thus |diocai|)- For example, a survey with ±15 deg Galac¬ 
tic cut with 30 million sources can detect the kinematic 
dipole at 5-(T, but needs to have a median redshift of at 
least Znied = 0.75 in order for this not to be excessively 
biased due to local structures. 

On the whole, Eig. [2] indicates that the convincing 
detection of the kinematic dipole expected given the 
CMB measurements is entirely within reach of future sur¬ 
veys, as long as those surveys have good coverage over 
both hemispheres and, given the source density, are deep 
enough not to be biased by the local-structure dipole. All 
requirements can be straightforward ly satisfied by sur - 
veys like some combination of LSST (jlvezic et al.ll2008l) , 
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Fig. 2. — Left panel: Statistical error in the dipole amplitude, marginalized over direction and other multipoles that are coupled to the 
dipole, as a function of the number of galaxies in a survey. The top horizontal dashed line shows the amplitude expected based on the CMB 
dipole measurements (A = 0.0028), and is the fiducial value in this work. The two dashed horizontal lines show the 3 (j and 5 (j detection 
of dipole with the fiducial amplitude. Right panel: defined in Eq. m, corresponding to the bias from the local-structure dipole, as a 

function of the number of objects A^gai- For a fixed amplitude of diocal Ax^ still depends on the direction of this vector; here we show 
the value averaged over all directions of diocal- Solid lines show cases when the median galaxy redshift is ^med = l-O? while dashed lines 
are for = 0.75. The legend colors are the same in the left panel. 


Euclid (|Laureiis et alJ[20Tlll and DESI (|Levi et alJl2013[l 
and, especially, by deep, all-sky surv eys with good red - 
shift information such as SPHEREX (|Dore et alJl2014ll . 

Einally, we have also calculated the statistical error in 
the direction of the kinematic dipole, based on the fidu¬ 
cial amplitude we had adopted. The direction’s error is 
generally rather large, e.g. an area of about ~ 10 deg in 
radius for Ng^\ = 10® and the Galactic ±15 deg cut. Nev¬ 
ertheless, a combination of the kinematic dipole’s ampli¬ 
tude and direction that roughly match the CMB dipole 
would present a convincing confirmation of the standard 
assumption. One could further carry out detailed fore¬ 
casts of what various findings could rule out the null 
hypothesis; we leave that for future work. 

4. CONCLUSIONS 

We have studied the prospects for measuring the kine¬ 
matic dipole - our motion through the LSS rest frame - 
as revealed by the relativistic aberration of tracers of the 
large-scale structure. The standard theory predicts that 
the kinematic dipole should agree with the CMB dipole, 
but this expectation could be violated due to a number 
of reasons. Therefore, verifying the standard expectation 
is an important null test in cosmology. The challenge 
comes from the fact that the dipole amplitude is small 
{A ^ 0.003), and easily contaminated by the intrinsic 
clustering of galaxies (the “local-structure dipole”). 

A successful measurement of the kinematic dipole 
therefore has two qualitatively different requirements: 
the survey should cover most of the sky and have enough 


objects to have sufficient signal-to-noise to detect the 
aberration signature of the dipole, but it should also be 
deep enough, so that the local-structure dipole contam¬ 
ination is sufficiently small. The two requirements are 
displayed in the two panels of Eig. [2] respectively. Eor 
a 5-cr detection, a survey covering > 75% of the sky 
in both hemispheres (our “Galactic ±15deg cut” case), 
with iVgai ± 30 million galaxies, is required. Eor a negli¬ 
gible bias, this same survey should have median redshift 
greater than about 0.75 or higher, with increasing depth 
requirements as Vgai increases. 

Eortunately these requirements can be satisfied by 
upcoming surveys, including DESI, Euclid, and LSST 
if they are combined, and potentially with SPHEREX 
alone. Even current all-sky surveys such as WISE (Wide- 
held Infrared Survey Explorer, iWright et al.l (|2010ll l 
are not out of the question, provided a sufficiently 
deep sample can be selected photometrically; current 
WISE samples hav e typical galaxy redshifts Zmed — 0.2 
(|Bilicki et al.ll20l4 l and are not yet deep enough to mea¬ 
sure the kinematic dipole. 
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